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Abstract 

The temperature dependence of the reflectance spectra of magnesium oxide (MgO) and fer- 
ropericlase (Mgi-^Fe^O, for x = 0.06 and x = 0.27) have been measured over a wide frequency 
range (~ 50 to 32 000 cm -1 ) at 295 and 6 K. The complex dielectric function has been determined 
from a Kramers-Kronig analysis of the reflectance. The spectra of the doped materials resembles 
pure MgO in the infrared region, but with much broader resonances. We use a shell model to 
calculate the dielectric function of ferropericlase, including both anharmonic phonon-phonon in- 
teractions and disorder scattering. These data are relevant to understanding the heat conductivity 
of ferropericlase in the earth's lower mantle. 

PACS numbers: 63.20.Kr, 63.20.Mt, 74.25.Kc, 78.30.-j 
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I. INTRODUCTION 



Ferropericlase, (Mgi-^Fe^O, with x = 0.10-0.15), is thought to be one of the major con- 
stituents of the earth's lower mantle (660-2900 km depth).- The name 'magnesiowiistite' 
is also used, but properly refers to the doping region x close to the wustite (x=l), rather 
than the periclase (x=0) limit. Transport properties of (Mg,Fe)0 are therefore important 
in modeling the Earth's thermal state and evolution, where both conduction and convection 
are operative.- 1 ^ The two heat carriers in conduction for an insulating mineral like ferroper- 
iclase are phonons and photons. Phonons (extended or localized) are distributed in the far 
and mid-infrared frequency range. They can be scattered by various defects (e.g. impuri- 
ties, grain boundaries, . . . ) and by the intrinsic anharmonic phonon-phonon interactions. 
Photons are described by Planck's black-body radiation formula, reaching energies ~ 1 eV 
at 2000-3000 K. The photon-matter interaction is more complex. Far- and mid-infrared 
photons usually couple with infrared-active vibrations. Photons of higher frequency can 
induce pure electronic transitions or vibronic transitions, depending on the details of the 
system's electronic states and adiabatic potential-energy surfaces. Goncharov et ai— mea- 
sured the optical absorption spectra of Mgi^Fe^O (with x=0.06, 0.15, and 0.25) across 
the high-spin/low-spin transition,-^ which occurs over a pressure range of 40-60 GPa at 
room temperature. Their results indicate that low-spin (Mg,Fe)0 will exhibit lower (rather 
than higher-) radiative heat conductivity than high-spin phase due to the red-shift of the 
charge-transfer edge. The origin of this spin transition and its influence on the radiative 
heat conductivity of ferropericlase are further investigated recently. 9 A complete picture of 
the thermal conductivity must include contributions from both phonons and photons. 

As a solid solution, ferropericlase has a vibrational frequency spectrum similar to that of 
pure MgO. However, with strong disorder scattering of propagating vibrational states, the 
harmonic eigenstates of the disordered crystal do not necessarily have a well defined wave 
number, and may not propagate ballistically. In addition, the anharmonic phonon-phonon 
interactions causes a shoulder at ~ 640 cm" 1 in the infrared (IR) reflectance spectrum of pure 
MgO. 10 Thus anharmonicity should also be included in analyzing the infrared reflectance of 
ferropericlase. 

In this paper we report the temperature-dependent infrared reflectance measurements of 
magnesium oxide and ferropericlase for several Fe concentrations at ambient pressure. We 
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construct a model in which anharmonic phonon-phonon interactions and disorder scattering 
are treated separately. Their effects are then combined for comparison with the experimental 
data. 



II. EXPERIMENTAL MEASUREMENTS 

The samples we examined are homogeneously doped single crystals, in which 
Fe 3+ / Yl F e ~ 0.02 for the 6% sample and 0.05 for the 27% sample. However, in our analysis 
the influence of Fe 3+ and magnesium vacancies is ignored. A detailed description on the 
samples' synthesis, crystallography and elastic properties is in Ref. ll|. The samples are 
rectangular slabs with typical dimensions of 1 mm x 2 mm, with a thickness of ~ 0.3 mm. In 
order to reduce interference effects due to reflections from the back surface, the samples have 
been wedged. However, due to the thin nature of the samples, the largest wedge that could 
be introduced was ~ 15°. The reflectance spectra has been measured at a near- normal angle 
of incidence at 295 and 6 K over a wide frequency range from « 50 to about 32 000 cm -1 
on Bruker IFS 66v/S and 113v spectrometers using an in-situ evaporation technique. 12 The 
measured reflectance at 295 and 6 K of pure MgO, and Mgi-^Fe^O, for x = 0.06 and 
x = 0.27 are shown in Figs. [H(a), (b) and (c), respectively. Although wedging the samples 
has been very effective at reducing interference effects, weak fringes may still be detected 
at low temperature below about 150 cm -1 . The complex dielectric function e = e± + ie^ 
has been determined from a Kramers-Kronig analysis of the reflectance, where extrapola- 
tions are supplied for to — > 0, oo. At low frequency, an insulating response is assumed and 
R(u -> 0) ~ 0.27, 0.28 and 0.31 for MgO, and the 6% and 27% Fe-doped materials, re- 
spectively. Above the highest measured frequency the reflectance has been assumed to be 
constant to approximately 75 000 cm -1 , above which a free-electron approximation has been 
assumed (R oc 1/cu 4 ). The imaginary part of the resulting dielectric function at 6 and 295 K 
of pure MgO, and Mg^^-Fe^O, for x = 0.06 and x = 0.27, are shown in Figs. [Tj(d), (e) and 
(f), respectively. The imaginary part of the dielectric function contains most of the physical 
information, and is the focus of our theoretical analysis. The optical features have been fit 
to a classical oscillator model using the complex dielectric function 

e(u J )=e oc + J2- - 9 , (1) 
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FIG. 1: The measured reflectance R(to); (a) pure MgO, and Mgi_ x .Fe x O for (b) 6% and (c) 27% 
Fe-doping. The corresponding imaginary part of the dielectric functions e2(w); (d) pure MgO, and 
Mgi_ x Fe x O for (e) 6% and (f) 27% Fe-doping. The solid line corresponds to data measured at 
6 K, dashed line corresponds to data at 295 K. 

where eoo is a high-frequency contribution, and lotoj, and u Pt j are the frequency, full 
width and effective plasma frequency of the jth vibration. The results of non-linear least- 
squares fits to the reflectance and €2(00) are shown in Table [H In addition to the strong 
feature in €2(0;) seen at about 400 cm -1 , other features at ~ 520 and ~ 640 cm -1 are also 
clearly visible in 62 (w) shown in Fig. [TJ however, these features are very weak and as a result 
the the strengths and widths of these modes are difficult to determine reliably. 



TABLE I: A comparison of the fitted values of the static and high-frequency contributions to the 
real part of the dielectric function at room temperature, as well as the fitted frequency, full width 
and effective plasma frequency (wto> 27 and u p , respectively) of the feature associated with the 
strong TO mode in MgO, and the 6% and 27% Fe-doped materials at 295 and 6 K. The units of 
lotOi 27 and top are in cm -1 . The strength of the TO mode is also expressed as a dimensionless 
oscillator strength S = uj^/uj^q. 
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"Values at 295 K, the estimated uncertainty is about ±0.1. 
b The uncertainty in wto is ±0.1 cm -1 . 

c The uncertainties for 27 are ±0.1 cm -1 in the pure material, and ±0.5 cm -1 in the Fe-doped materials. 
d The uncertainty in u> p is ±20 cm -1 . 



III. COMPUTATIONAL METHODS 



A. General Scheme 



Infrared dielectric properties of ionic crystals are contained in the linear response function 
£ap(v) = e aj a(oo)+47rx a/3 (u;).— Considering only the first-order moment of the electric dipole, 
the dielectric susceptibility of a crystal can be related to its displacement-displacement 
retarded Green's function by: 

XafiM = "ttfV f^^^([D a (t),D p (lf)])e^U(t-lf) 
NV C > /_ 00 ih 

" 9{t ~ l ' ([u 7 (ls; t),u s (l's'; t')])e i " ( '-'' ) d{t - t') 



c Is-y l's'5 J ~°° 

\f E E Z a7 (ls)Zp S (l's')G yS (ls, l's'- oo), (2) 



NV C 

Is-y l's'5 

where D a (t) = Z a p(ls)up(ls; t) is the a component of the first order electric dipole of 

lsf3 

the whole crystal, Z a p(ls) is the Born effective charge tensor of the atom s at site I, and 
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up(ls; t) is the atom's displacement at time t. The volume of a single cell is V c , and iV is the 
number of the cells in the whole crystal. The Green's function G a p{ls, I 1 V; t — If) is defined 
as: 

G a p(lsJ's';t-t') = 9 JL^l([ Ua (l S]t ),up(l's'-,t')}), (3) 

which can be evaluated from its equation of motion. 14 For a harmonic crystal, the vibrational 
Hamiltonian is quadratic and can be solved exactly. We denote the eigenvectors of a pure 
crystal as -^e Q ,(s|qj)e iq ' R ^ s \ the corresponding eigenvalues as Uqj, those of a disordered 
crystal as e a (s\j) and Uj, the Green's function of the pure as g, the disordered as G°. Then 

e a (s|qj-)e^( S '|q7-)e^W s )- R ^')) 

Q a B(ls, i s :u>) = > = , (4) 

V iVv/M(s)M( S ')(^ 2 - + ^) 

G£ fl (Zs,ZV;u;) = V ; V ' ; ^ — — , 5 

P Y ^M{ls)M{l>s'){uj 2 - J] + t2u V ) 

where the mass of the atom s is denoted as M(s) in the pure crystal, M(ls) in the disordered 
crystal, with the extra label I to specify its site, rj is an infinitesimal number ensuring 
causality. 

Anharmonic interaction will couple these modes and make exact solution impossible. The 
standard treatment of this many-body effect uses the Dyson equation to define a self-energy 
for each mode. We can either choose e a (s\j) as the unperturbed states, then the only inter- 
action will be anharmonicity, or choose ^^e Q ,(s|qj)e* q R ' ( ' /s - ) as the basis and treat disorder 
as an extra perturbation. The first approach has been used by one of the authors (PBA) to 
study the anharmonic decay of vibrational states in amorphous silicon.— In this paper we 
use a hybrid approach. We write the dielectric function of a disordered anharmonic crystal 
in the perfect crystal harmonic basis as 



e a(3 (uj) = e a(3 (oo) +47rx Q/3 M 



,e 7 (s|0j)v- 7 f^Wj) 



T() ^ VW*)^ P WW) 



= e Q/3 (oo) + y c 2^ M0j) 2 _ U 2 + 2 W (oj)(A(Oj,^) - ir(Oj»}' (6) 

where ou(Oj) = ct>TO,j is the frequency at q = of the jth TO branch. The terms A(Oj,u) 
and T(Oj,u) correspond to the real and imaginary part of the mode's self-energy S(Oj, uj). 
Then we split this self-energy into two parts: £ = S anharm + S disorder . Each piece is cal- 
culated independently. This is equivalent to omitting all the diagrams where the disorder 



scattering vertex appears inside an anharmonic interaction loop. The anharmonic interac- 
tion of ferropericlase is assumed to be the same as that of pure MgO, i.e. the influence of 
disorder on anharmonic coefficients is totally ignored. Disorder is treated by exact diago- 
nalization without anharmonicity, it is then converted to a self-energy of the TO mode in 
the perfect crystal harmonic basis. These approximations are tested by comparing with the 
experimental results. 

B. Shell Model 

The scheme described above is general. It does not depend on which microscopic model 
is chosen to get harmonic phonons, disorder scattering strength, and higher-order force con- 
stants. Here we use an anharmonic shell model, with shell parameters fitted to experiments. 
The harmonic phonon properties in this paper are calculated with the general utility lattice 
program (GULP) code. 16 

Two sets of shell parameters' 6 - 1 ^ are used for MgO: S-I and B, and one for FeO: S-II. 
S-I and S-II are rigid shell models in which O 2- has the same set of parameters, thus they 
can be conveniently used to simulate ferropericlase. B is an isotropic breathing shell model 
which gives better fit to the experimental data. However, it can not be directly used for 
ferropericlase. For FeO the elastic constants C12 >C44, while the isotropic breathing shell 
model is only suitable for cases where C12 <C44.— We treat B as a reference to check 
our anharmonic calculations based on S-I. All the model parameters are listed in Table HT1 
Table HTT1 contains the calculated physical properties and corresponding experimental values. 
Phonon dispersion curves for the pure crystals of MgO and FeO are shown in Fig. [2j 

C. Anharmonicity 

A complete calculation of anharmonicity is tedious, even for a pure crystal.™ Thus we 
ignore the less important terms and focus on the dominant one. From Eq. ([6]) it is clear 
that since |E] = |A — iT\ is small compared to c^to, the real part of the self-energy A has 
negligible influence on t2{uj), except to shift its resonant frequency. The shell models we use 
are fitted to the experimental data measured at room temperature. The anharmonic shift is 
small, compared with the shift caused by impurity scattering. Thus, we ignore it completely 



7 



TABLE II: Shell model parameters used in the calculation . > The short-range repulsive potential 
is assumed to be a two-body Buckingham type: for S-I and S-II, V(r) = Ae~x.p(—r / p) — C/r 6 ; for 
B, V(r) = j4exp(— (r — ro)/p) — C/r e . The parameter k represents the spring constant between 
core and shell. Rows in which atomic symbols have a star (*) are for the B model. The label 
'shell' denotes a potential that acts on the central position of the shell, while 'bshell' denotes an 
interaction that acts on the radius of the shell which was fixed at 1.2 A. An extra parameter in B 
model is kesM = 351.439 eVA -2 . The equilibrium shell radius tq is 1.1315 A after optimization. 
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TABLE III: Physical properties of pure MgO and FeO, compared with shell model results. 
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FIG. 2: Phonon dispersions of the pure crystals, (a) MgO, solid line corresponds to the rigid-shell 



model S-I, dashed 
taken from Ref. 



ine to the isotropic breathing-shell model B, dots are the experimental data 



18 (b) FeO, solid 



experimental data taken from Ref. 
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and only consider the imaginary part of the self energy r an h ar m(0j, uj). To the lowest order 
r a nharm(0j, oj) can be written as^ 
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{(ni + n 2 + l)[S((Vi + u 2 — co) — 8{ux + uj 2 + uj)] 



+ (n 2 - ni)[8(u2 —uji + lo) — 8(u\ - uo 2 + to)]}, 



(7) 



where n\ = n(qji) is the Bose-Einsein population factor of the mode, and U\ = u(qji) is 
the corresponding frequency. The anharmonic coefficient V(°. ^ is 

'° q ,~ q )4£££*< 



V 



J Jl 32 



.„ / / J e Q (si |0j)e^(s 2 |qji)e 7 (s 3 | - qj 2 ) 
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i Gxp \iq- \H(l 2 s 2 ) — Hfhs^)}} . (8) 

The third-order force constants $q,/3 7 ( a ° ' s 2 2 ^ ) are large only for the nearest neighbors. Sym- 
metry will restrict most of them to be zero, and among those nonzero terms only two are 
independent. The general formula for third order force constants is^ 

'o o r 
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BR a RfjR^ + C^RaSp^ + R^8 a -y + R^8 a fj) 



B R 3 R* 



C 
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(9) 



where R is the lattice distance between the ion (°) and ([,) , and R a is its projection along a 
direction. The term <p(r) is the two-body pair potential, and <f)', <$"... are derivatives with 
respect to r. Following E. R. Cowley," we compute $ a /3 7 (° by direct differentiation over 
the nearest-neighbour short-range potentials and Coulomb potentials. For the rigid-shell 
model S-I, (f>(r) = Aexp(—r/p) — — . For the breathing-shell model B, <j)(r) = v4exp(— (r — 
r o)/p) — ~- To be more specific, if we take a Mg 2+ as the origin and denote it as 1, its 
nearest neighbor 2 ~ along the [100] direction as 2, then from symmetry we can determine 
$ xxx (112) = $ xxx (121) = -$ xxx (122) = $ yyy (112) = $ xyy (112) = $ xzz (112) 
Putting in numbers from Table U we obtain $ xxx (112) = -25.34 eV/A 3 , $ xyy (112) = 
-1.79 eV/A 3 for the S-I model, and $ xxx (112) = -24.0 eV/A 3 , $ xyy (112) = -1.78 eV/A 3 for 
the B model. If we do not include the Coulomb interaction, these values will be $ xxx (112) = 
-42.70 eV/A 3 $ xyy (112) = 6.88 eV/A 3 for the S-I model, and $ xxx (112) = -41.57 eV/A 3 , 
$ xyy (112) = 7 eV/A 3 for the B model. Although $ xyy (112) is small compared with $ xxx (112), 
it can still have non-negligible influence on the amplitude of the T(0j,u) near 640 cm -1 . 
Other parameters (Born effective charge tensors, harmonic eigenvectors) are obtained from 
GULP. The integration over q-space is done with the tetrahedron method, using 1/48 of the 
Brillouin zone, and averaging over x, y, and z polarizations. We use 3345 q-points, equivalent 
to 160560 q-points in the whole Brillouin zone. 



D. Disorder Scattering 

The self-energy of a vibrational mode caused by disorder scattering is defined 
statistically,^^ 

«G°»=g + g£«G°», 

where ((G )) denotes the Green's function averaged over different impurity distributions. 
We slightly modify this definition by including the Born effective charge. From Eqs. (j4j) and 
we define the following equation 

«E E Zay(ls)G° yS (ls, IV; u)Z S p(l's'))) = E E ZMMl*, l'*t\ ^)Z S p(s'), (10) 

7 8 7 5 

where 

^ e 7 ( g |qJ■)e|(y|qJ^)e ^ q•( R ^- R (^ , )) 

q~s[ls,l s :u>) = > = . (11) 

V N^M{s)M{s>){^ - ojIj - z2^S disordcr (qj- v)) 
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The self-energy defined in this way guarantees that the dielectric susceptibilities calcu- 
lated from G° and g are the same. Summing over all sites of the crystal leaves only TO 
modes on the right hand side of Eq. (1101) . Thus, once we get the averaged dielectric sus- 
ceptibility ((Xa/?)) from the exact eigenvectors of the disordered crystal, we can extract the 
self-energy of its TO phonon. 

We expand an orthogonal 8-atom MgO unit cell in each direction by 5 times, which gives 
a 5 x 5 x 5 super-cell containing 1000 atoms. Then we randomly replace the corresponding 
number of Mg 2+ by Fe 2+ . The shell parameters of Mg 2+ are from S-I model, those of Fe 2+ 
are from S-II model, those of O 2- are the same in both models. From Eqs. (j2J) and (jHJ), for 
each configuration we have a harmonic susceptibility 

**H = nv c ^ u 2 - u 2 - i2un • (12) 

3=1 J 

We can choose a small value for rj and evaluate Eq. ffT2]) directly (Lorentzian broadening). 
However, insofar as rj is finite, it is equivalent to have each mode j in Eq. (1121) an imaginary 
self-energy (life time) linear in frequency u. The self-energy of the TO phonon ^disorder 
extracted from this approach will depend on frequency linearly. Replacing the factor 2u by 
2ujj won't help either, as each mode j now has a life time independent of frequency, and 
^disorder will be a constant depending on 77 when u — > 0. To avoid such artifacts we use 

i I Z7T 

+ — [6(u - Uj) + 6(u + Uj)] 



u^j — u 2 — i2ur] u 2 — u 2 2u 

to separate the real (xi) and imaginary part (^2) of the dielectric susceptibility. Then 
we divide the vibrational spectrum into equally sized bins (1 cm" 1 ) and compute \2 as a 
histogram. The real part xi is obtained from \2 from the Kramers- Kronig relation. Many 
such super-cells are built and their eoo and x calculated. We find that 10 configurations 
are sufficient to give a well converged average. The final and x are assumed to be the 
averaged values of all configurations. To remove the unphysical spikes caused by the finite 
size of our super-cells, while keeping the main features unchanged, we further smooth the 
dielectric susceptibility by averaging over adjacent bins iteratively, 

xT x iJ) = 1 1x50' - 1) + 4x50') + xl{j + 1)] • (is) 
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FIG. 3: Computed anharmonic properties compared with experimental data for pure MgO. (a) The 
imaginary part of dielectric function at 295 K; the experimental data are the same as those in 



Fig. 
Ref. 



rjd). (b) The imaginary part of self energy at 295 K; the experimental data are digitized from 
25l . which are fit to infared spectra based on a semi-quantum dielectric function model. 



In this way we successfully simulate the dielectric function of a 'real' crystal (real in the sense 
that except for finite size, disorder scattering is treated without any further approximations). 
These results, together with anharmonicity, are summarized in the next section. 



IV. COMPARISONS AND DISCUSSION 



The anharmonic effects in pure MgO are shown in Fig. [3j The computational results and 
experimental values are quite close, especially near 640 cm -1 which corresponds to TO+TA 
combination mode. It is not surprising, since the procedure we followed was originally 
developed and worked well for alkali-halide salts, which are similar to MgO. Below 800 cm"" 1 , 
the rigid shell model S-I and breathing shell model B give almost identical self-energies. The 
discrepancy in the high-frequency range indicates that the dispersion relations from empirical 
models are less accurate for high-frequency optical branches. The width of TO mode at the 
reststrahlen frequency c^to is less accurate, as r anharm (0j, uj to ) is intrinsically small and 
higher-order anharmonic effects become important.— 

Figures H] and [5] show how anharmonicity and disorder scattering influence the dielectric 
function. For the 6% sample it is clear that the shoulder near 640 cm -1 is caused by 
anharmonicity, while the shoulder at about 520 cm -1 is due to disorder scattering. Disorder 
scattering becomes stronger for the 27% sample and seems contributes to all the shoulders. 
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FIG. 4: The anharmonic and disordering scattering effects in Mgi_ x Fe x O for the 6% Fe-doped 
sample, (a) Imaginary part of the dielectric function. The labels '6 K Calc' and '295 K Calc' 
denote the calculated curves, including both disorder scattering and anharmonic interactions at 
the corresponding temperature. Experimental data are the same as those in Fig. Q^e). The label 
'No anharm' denotes the dielectric function calculated from disorder scattering only, (b) Imaginary 
part of self energy. The labels '6 K' and '295 K' denote the self-energies caused by anharmonic 
interaction at the corresponding temperature; 'disorder' denotes the self-energy due to disorder 
scattering, which is computed by histogram method where the bin size equals 1 cm -1 , then itera- 
tively averaged 30 times. The total self-energies are the sum of these two pieces, and are used in 
calculating the '6 K Calc' and '295 K Calc' dielectric functions shown in (a). 

The shoulder caused by anharmonicity corresponds to a peak in the two-phonon DOS, while 
shoulders caused by disorder scattering are related to peaks in the one-phonon DOS. 

Figure [6] contains the reflectance computed from the dielectric functions at 295 and 6 K 
shown in Figs. H] and [51 As in the case of pure MgO, the agreement between theory and 
experiment is better in the region where the self-energy caused by lowest-order pertubation 
is large. Near the reststrahlen frequency cjto, the self-energy is smaller, and RiiS) is more 
sensitive to details. Our model underestimates the broadening of the resonance, but correctly 
identifies the sources of broadening. 

It is of interest to determine whether the disorder scattering is mainly due to the differ- 
ences in mass or in the inter-atomic potential. Thus we repeat the above procedure with 
a model which only contains mass disorder, i.e. Fe is treated as an isotope of Mg, its shell 
parameters are the same as Mg 2+ in S-I model. It turns out the most significant factor 
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FIG. 5: The anharmonic and disordering scattering effects in Mgi_ x Fe x O for the 27% Fe-doped 
sample, (a) Imaginary part of the dielectric function; (b) Imaginary part of self energy. The 
computation procedure is the same as for the 6% Fe doping. 




FIG. 6: The calculated infrared reflectance, compared with the experimental data (same as in 
Fig. mt>) and (c)) for Mg^Fe^O. (a) 6% Fe doping; (b) 27% Fe doping. 

is : For the isotope model (S-I) is the same as pure MgO (2.94), for S-II model 6^ 
increases to 3.05 for 6% Fe and 3.47 for 27%, in reasonable agreement with the results shown 
in Table [H Thus the LO frequency predicted from the isotope model is larger than the ex- 
perimental value. The differences in the inter-atomic potentials change the relative strength 
of the self-energy, but in both cases the self-energy spectra carry features of the one phonon 
DOS of pure MgO. 

In addition to phonons, electronic transitions may also influence the infrared dielectric 
properties of ferropericlase. Wong— measured the far-infrared absorption spectra of iron- 
doped MgO. A line at 105 cm -1 was observed with a peak absorption coefficient of 1.5 cm -1 
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and a width of ~ 9 cm -1 at 20 K in a sample with 0.2% Fe. This feature is attributed to the 
transition r 5g — ► r 3g , r 4g of MgO: Fe 2+ at cubic sites. If we assume the absorption coefficient 
is proportional to the impurity concentration, then we can estimate the corresponding e2 



105 cm -1 ), a(u>) is the absorption coefficient at frequency u (in units of cm -1 ). The value of 
€2 is about 0.22 for 6% Fe concentration, 0.98 for 27%. As the iron concentration increases, 
the electronic transitions of Fe 2+ should show greater influence on the far-infrared spectra 
of ferropericlase. In our measurement the spectra below 200 cm -1 are complicated due to 
the presence of fringes, consequently we can not confirm this tendency. Henning et al. 27 
measured the infrared reflectance of FezMg^zO for x = 1.0, 0.9, 0.8, 0.7, 0.5 and 0.4 at 
room temperature. The €2 curves reported in their paper do not show a monotonic rise in 
the far-infrared region as the iron concentration x increases from 0.4 to 1.0, while they are 
all in the range of 6-10 near 100 cm -1 . It is difficult to explain such large €2 with lattice 
vibrations alone, and the accuracy of these data has been questioned.— Further experiments 
are needed to clarify this issue. 

V. CONCLUSIONS 

The infrared reflectance spectra of magnesium oxide and ferropericlase has been measured 
at 295 and 6 K. It is found that increases as Fe concentration increases, while the width 
of the TO modes remains the same in the doped materials. We construct a theoretical model 
which includes both disorder scattering and anharmonic phonon-phonon interactions. The 
model shows fairly good agreement with the experiment in the regions where the lowest- 
order perturbation is relatively large. Near the resonance, theory and experiment both have 
smaller self-energies, which makes the reflectance quite sensitive to the details. We do not 
know whether the disagreements with experiment in the region are caused by neglect of 
higher order corrections, or by inaccuracy of the underlying model. However, the model 
identifies the global features reasonably well, and may provide a good basis for the study of 
phonon decay needed for a theory of heat conductivity. 




where n is the refractive index (for pure MgO, n ~ 3.2 at 
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